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An analytical study is taken to investigate the relationship 
between material softening and structural softening through 
the use of a model problem in one dimension. With general 
nonlinear assumptions on the constitutive relations, it turns 
out that the governing equations can be viewed as a system 
of parametric equations, which couple the size effect and the 
nonlinear effect. Compared with the bilinear assumptions in 
previous literature, we find that the nonlinear assumptions 
herein capture more details in the post-peak structural re- 
sponse. After doing standard mathematical analysis to the 
nonlinear equations, we manage to derive necessary and suf- 
ficient conditions for the occurrence of four important post- 
peak cases, which are often observed in experiments. In par- 
ticular, our analysis reveals that the mechanism of the snap- 
through phenomenon is due to the convexity change of the 
constitutive curve of the softening part. Mathematical exam- 
ples are also given to illustrate the proposed procedures. 



1 Introduction 

Strain-softening, i.e., the decrease of stress with the in- 
crease of strain, is such a common phenomenon that has been 
recorded for a variety of materials, like concrete, rocks, ce- 
ramics, metals, etc. Bazant et al., [1] gave a comprehensive 
review of this phenomenon and analyzed its mechanism from 
a continuum point of view. Moreover, it is well-known that 
strain softening is always accompanied by highly localized 



deformations of the specimen ( [2, 3|). Due to the impor- 
tance of softening phenomenon in structural safety assess- 
ment, many efforts have been made in the past decades to 
investigate strain-softening with localization experimentally, 
numerically, and analytically, as reviewed by flUE] • 

Snap-back may be one of the most interesting and per- 
haps most common structural instability phenomena ob- 
served in experiments. It shows that the load-displacement 
curve displays a positive slope after attaining the peak load, 
de Borst |6| demonstrated the possibility of snap-back be- 
havior on structural level by means of two concrete struc- 
tures: a reinforced concrete and an unreinforced specimen. 
In order to simulate the highly localized failure mode in a 
strain-softening solid, a modified arc-length control method 
was used in that paper. Later, Rots and de Borst |7| did a 
tensile test on concrete specimens and analyzed it by using 
the finite element method, with a particular attention on the 
snap-back behavior. He et al., O studied the class II be- 
havior (snap-back) of rock with a spring model, which was 
characterized by non-uniform failure. Unloading-reloading 
tests were also conducted in the post failure region in that 
paper. One of their results is that, if inelastic strain increases 
slower than the elastic strain decreases, rock shows class II 
behavior. 

Jansen et al. J9| did an experiment on concrete cylin- 
ders by using the feedback-control method. From two test 
series, the stress-displacement behavior for different height- 
diameter ratios with normal strength and high strength were 



obtained. They found that the pre-peak segment of the stress- 
displacement curves agrees well with the pre-peak part of 
the stress-strain curves, while the post-peak segment shows a 
strong dependence on the geometric size, namely the radius- 
length ratio. More specifically, the longer the specimen is, 
the steeper the post-peak segment of the stress-displacement 
curves becomes. The feedback-control method was also used 
in Subramaniam et al. [ 1 1 to test concrete in torsion, and 
snap-back was also found in the experiment. 

Some analytical studies were also taken to investi- 
gate softening with localization. With the use of a one- 
dimensional model, Schreyer and Chen ifTTI analyzed the 
snap-back phenomenon and found the important size effect 
on the instability. Due to the simplicity of bilinear assump- 
tions on the constitutive relations, further features like snap- 
through were lost in the result, although in some experiments 
this feature was observed (see van Vilet and van Mier 1 12 1). 
The same constitutive relations were also assumed in Chen 
et al. lfT3ll to analyze the stability in some hierarchical struc- 
tures. In a more complex setting with certain nonlinear as- 
sumptions on the constitutive relations, Sundara Raja Iyen- 
gar et al. 1 14 1 took an analytical study. By using the fictitious 
crack model (FCM) developed by Hillerborg, they found the 
effect of the softening exponent n on the size effect and snap- 
back behavior of beams, while the stress-displacement rela- 
tion was assumed as a general power law function. Dai et 
al. Ifl5ll constructed the analytical solutions for localizations 
in a hyperelastic slender cylinder. With the use of coupled 
series-asymptotic expansions approach and phase plane anal- 
ysis, they solved the partial differential equations and found 
that the width of the localization zone depends on the mate- 
rial parameters in the post-peak region. Further, they showed 
that there is a snap-back phenomenon when the radius-length 
ratio is relatively small, which agrees well with experimen- 
tal observations. Dai et al. Q~6| showed a similar result for 
hyperelastic shape memory alloys. Gradient theory may be 
another powerful tool in dealing with localization of defor- 
mation (see Triantafllidis and Aifantis |17|). For example, 
Triantafllidis and Bardenhagen lfT8l investigated the issues 
of instability and imperfection sensitivity of the solutions of 
a boundary value problem in one dimension. Their results 
also revealed some important size effect. 

To the authors' knowledge, however, there is not any an- 
alytical study with general nonlinear constitutive relations in 
the open literature which explores the role played by the con- 
vexity of the constitutive curve of the softening part and the 
coupling effect between this convexity and the size. Also, 
both snap-back and snap-through were observed in some 
experiments, but no analytical results are available for ex- 
plaining the transition from snap-back to snap-through. We 
shall explore these aspects in this paper. To gain insight into 
the post-peak response, we study the same one-dimensional 
structure as considered in lllllll91l20ll2D . The difference is 
that here we use general nonlinear constitutive relations, in- 
stead of the bilinear ones used in these papers. First, we set 
up the stress-strain equations for the structure in the post- 
peak region, which are nonlinear as compared with the bi- 
linear case. After some analysis, we derive the mathemati- 



cal conditions for the occurrence of several important curves 
as frequently observed in experiments, including the snap- 
through (which cannot be captured by the bilinear assump- 
tions). Finally, an example is given to illustrate these cases, 
and the post-peak curves are consistent with our theoretical 
predictions. 



2 Model Problem 

To simulate post-peak experiments, we consider a struc- 
ture with a serial arrangement of intact elastic and strain- 
softening zones. This model was used by several researchers, 
such as 1 19 20, 21] in the early years. In | fTT) . it was in- 
troduced to analyze strain-softening with bilinear assump- 
tions on the constitutive relations. As shown in Figure 1(a), 
the structure is a bar of length L = a + b with a unit cross- 
sectional area. That is to say, it is composed of two segments 
(segment A with length a and segment B with length b). The 
two segments are usually described by similar constitutive 
equations, and the main difference is that the limit stress for 
B is slightly less than that of A. Therefore, if the stress on the 
structure is such that the strain in region B exceeds the value 
at the limit state, then softening will occur. It is assumed that 
softening occurs uniformly over a localized region B under 
quasi-static loading. 

In order to consider a general nonlinear case, the consti- 
tutive relations for the two regions are set as: the loading and 
unloading segments of region A are two arbitrary functions 
/n and fa respectively, while the loading and softening seg- 
ments of region B are two arbitrary functions fa and fa re- 
spectively. Moreover, we assume that the foregoing nonlin- 
ear functions are twice differentiable with f n > 0, f 12 > 0, 
f 2i > 0, f 22 < 0. The limit stress for region A, denoted by 
a a , is assumed to be slightly larger than that of region B, 
which is denoted by Go. The details are shown in Figure 
1(b) (where f\ is used to denote both the pre-peak and the 
post-peak segments, as region A only experiences loading or 
unloading ). As the post-peak curve of the structure is our 
main concern, in the following derivation, for simplicity, we 
use /i , fi to denote the post-peak curves of region A and 
region B respectively, unless otherwise specified. 

As to the post-peak response, for a strain softening ma- 
terial with a serial setting (cf. Figure 1(a)), region A is in 
an unloading process and region B experiences strain soften- 
ing. Given the values of strain in regions A and B, say e\ and 
e2, respectively, then the composite strain for the complete 
structure is given by 



aei+be 2 .. . 
e = = (1 —n)e\ +ne2, 



(1) 



where n — b/L. Since we consider it as a quasi-static prob- 
lem, the composite stress is then given by 



a = /i(ei) =h{ei)- 



(2) 



In fact, one can easily see that (1) and (2) are also true if 




obtain 



G10 G 



(b) 



Fig. 1 . (a) One-dimensional model problem; (b) Stress-strain rela- 
tions for A and B. 



/ 2 is used to denote both the pre-peak and post-peak seg- 
ments. Here, for the post-peak region, we consider only 
when a >a* (a* represents the lowest stress value at which 
the bar breaks), and denote <?n and e 2 \ the values such that 
/l (en) = / 2 (e 2 i) = O*. Then, for the post-peak region, we 
have e\ £ [en, do] and e 2 6 [e 2 o,e 2 i] (see Figure 1(b) for 
the definitions of e\o and e2o)- From equation (2), we get 
e 2 =f 2 l [fi(e l )] (orei =/f 1 [/ 2 (e 2 )]). Thus ( 1 ) and (2) can 



be transformed into the system 



(l-ra)ei 



nf2 l [fM)]. 



(3) 



& = (i-*)+»4 



fl(ei) 
/jfe) 



(4) 



If (1 — n)f 2 (e 2 ) + nfi (ei) ^ 0, we have 

do = /j(ei)/ 2 (e 2 ) 

de (l-n)/ 2 (e 2 )+n/|(ei) ' 



</ 2 a _ «[/i( e i)] 3 / 2 >2) + (l -n)^(ei)[/ 2 (e 2 )] 3 



</e 2 



[»^( ei ) + (l-")/2( e2 )] 3 

In order to analyze the sign of (5), we define 

g(ei,er,n) = ra/i(ei)H-(l-n)/ 2 (e 2 ), 



•»(«!, d) = 



/ 2 («2) 



G(e u e 2 ) = [f 1 (e l )] 2 f 2 (e 2 )-[f 2 (e 2 )} 2 f 1 (e l ) 



(5) 



(6) 



(7) 



(8) 



(9) 



The above three functions can be viewed as functions of ei- 
ther e\ or e 2 by the relations between them as shown above. 
We note that »i(ei,e 2 ) depends on the slopes (the first-order 
derivatives) of the constitutive curves, G(ei,e 2 ) depends on 
the convexities (the second-order derivatives) of the consti- 
tutive curves and g(ei,e 2 ;n) depends on the size parame- 
ter n. We also point out that m(ei,e 2 ) = n is equivalent to 
g(ei,e 2 ;n) = 0. We shall see that n has an important influ- 
ence on the structural response. 



3 Post-peak Curves and Conditions 

Assuming that / 2 (e 2 o) = and f 2 (e 2 rj) < 0, that is to 
say e 2 o is a local maximum of f 2 . Then, for different f\, f 2 
and n the four cases shown in Figure 2 can arise. Next, we 
shall establish the conditions for each case. 

3.1 Case A: Stable Softening 

For the structure to be in stable softening (i.e., d<5Jde < 
0), from (5), it is easy to see the necessary and sufficient 
condition is 



which can be viewed as the parametric equations for the en- 
gineering stress-strain curve. We note that n is actually the 
width (scaled by L) of the localization zone in the reference 
configuration, as material points in region B are in the local- 
ization zone in the post-peak region. Obviously, system (3) 
couples the size effect and nonlinear effect. 

Now, we differentiate system (3) with respect to e\ to 



g{e u e 2 ;n) =«/ 1 (ei) + (l 
From which, we get 



■n)f 2 {e2) >0,forei e [«n,cio]. 

(10) 



— = n > no := max m(e\ , ei) - 

L «lG[en,«lo] 



max m(ei,e 2 ) . 

e 2 G[e20 ! e2l] 

(ID 



Case A 



Case B 




So, we get 



g(en,e 2 i;n) =n/ 1 (en) + (l-n)/ 2 (c2i) >0. (13) 



Thus, 



Since 



f 2 (en) 



> 



> -■ 



»o 



1 — n 1 — no ' 



(14) 



(15) 



Fig. 2. Four cases of the post-peak engineering stress-strain curves 



«<) 



3.2 Cases B and C: Snap-Through 

Now, we focus on the interval n £ (0,«o]- There are sev- 
eral possibilities, as shown in Figure 2. Before analyzing the 
remaining cases, we point out that the initial part (i.e., the 
part close to the peak) of the post-peak curve is in a state of 
stable softening for the conditions imposed on /] and f%. In 
fact, g(eio,e2o) = «/i(<?io) > , and at the peak point, we 
have da/de < 0. By continuity, there must be a part of the 
post-peak curve for e close to So (So = (1 — n)e\o + ne2o) 
in which da/de < 0. Also, at e\ = eio,<?2 = e 2 Q, we have 
d 2 a/de 2 = f 2 (e2o)/n 2 < 0. This would be useful for our 
later derivation. 

We see that each of Case B and Case C represents a 
snap- through case. Here, snap-through is defined to be the 
point at which the slope of the force-displacement curve be- 
comes infinite. As a result, when displacement (elongation) 
crosses this point, the force may experience a sudden drop. 
Firstly, let us consider the similarities between Case B and 
Case C. There are two turning points (the points at which 
da/de = oo) in both curves. From (5), it can be seen that this 
is equivalent to that the equation 

g(e u e2;n)=nf' l (ei) + {l-n)f2(e 2 ) = (12) 



has two roots, say e* n and e* 12 {e* n > e* l2 )- The following 
theorem provides a necessary and sufficient condition for the 
occurrence of the two turning points. 



we have 



mm -j 

1— »0 «ie[eii,«io]/i(ei) 



/i(en) *i €[«!!, *io] 7i(ei) 



(16) 



(17) 



Suppose that the minimum is attained at en (the corre- 
sponding e2 is given by / 2 - [/i(ei2)] = 622)- That is, 

/ 2 (e22)//l0i 2 ) = min {/ 2 (e 2 )//|(ei)}- If we view 

eie[eil,«10j 

f^\ e ^)l f\ ( e as a function of e 2 , then according to the La- 
grange Mean Value theorem, there exists an t, £ (e22,e2i) 
such that 



d r / 2 ( e z) i 
de 2 f[( ei ) 



\ei=\ 



/ 2 ( e 2l) _ / 2 («22) 

/l( f ll) /i'( e 12) 

e21 —^22 



>0. 



(18) 



As 



d r/ 2 M 1= [A(e 1 )} 2 f2(e2)-[f2(e2)] 2 f';(e i ) 

de2 [ fMY IfMW 



we have 



Theorem 3.1 . If two turning points occur, then the function 
G(e\,e2) must change sign at least once for e\ £ [en,eio]. 
On the other hand, if the sign ofG(e\, e 2 ) changes only once 
fore\ £ [cn,Cio]i then for any n £ [n\,no], two turning points 
occur, where n\ = m(e\\,e2\)- 

Proof. If two turnings occur, then the sign of the function 
g(e\,e2) changes twice (cf. Case B or Case C in Figure 2). 



[/1 [QYfi {%) - [fliWh (Q > 0, where C = /f 1 \fi{\)] , 

(20) 
which implies that G(^,4) > 0. Since G(eio,e2o) = 
[/i( e io)] 2 /2( e 2o) <0, the sign of G{e\,e 2 ) changes for e\ £ 
[eii,eio]. 

On the other hand, suppose that for e\ £ [en,eio], 
G{e\,ei) changes sign once. Now, we consider the function 
m{e\,e2) (cf., (8); we regard it as a function of ei). We now 



m 



While for Case C, we have 




81 = (l-nK 2 +n/ 2 - 1 [/ 1 (4 2 )]<So 



(24) 



In other words, in Case C the post-peak curve has entered 
the pre-peak region, while in Case B it has not. We find 
that, for given f\ and / 2 , there are some conditions for the 
occurrence of Case C. For simplicity, we assume that the sign 
of G(ei,e 2 ) changes once, say, 



e 2 o e 2n 621 

Fig. 3. Schematic illustration of m(ei,e 2 ) 



G(ei,e 2 ) <0, e 2 G [e20,e22), 

G(ei2,e22) =0, 

G{e\,e 2 ) >0,e 2 G (022,621] • 



(25) 



The following theorem provides a necessary and sufficient 
condition for Case C. 

show that for any n G [«i,«o], equation (12) has two roots. In Theorem 3.2. Under assumption (25) and n G [ni.no], a 
fact, it is easy to get necessary and sufficient condition/or the occurrence of Case 

C is 

dm _ -G(ei,e 2 ) 

de2 [/ 2 ( e2 )-/i( e i)] 2 /;(ei)' / 2 (e22)(42-«20)>/l(«12)(«12-«10)- (26) 



Thus, dm/de 2 also changes sign once. We also note that •'' *• •" 

«o is a maximum of m(ei,e 2 ), say, attained at e2«. Then 
dm/de 2 — at <?2 = <?2«- On the other hand, 



dm 

de~ 2 le2=e20 



-G(eio,e 2 o) 



[/ 2 (^o)-/i( e io)] 2 /i( e io) 



>0. 



(22) 



So, the curve m{e\,e 2 ) should have the characteristics shown 
in Figure 3. Thus, for any n G [nj ,no], n — m[e\ ,e 2 ) has two 
roots, which then implies that g(ei,e 2 ;n) has two zeros. This 
completes the proof of the second part of the theorem. 

If f\ is linear in the post-peak region, then G(ei,e 2 ) = 
L/i( e i)] 2 /2 ( ez )- Consequently, the sign of G(e\,e 2 ) depends 
on the convexity of f 2 . We have the following corollary: 

Corollary 3.1 . For f\ being linear in the post-peak region, 
if two turning points occur, then the convexity of f 2 must 
change at least once for e\ G [en,eio]. On the other hand, 
if the convexity of f 2 changes once, then for any n G [ni,no], 
two turning points occur in the post-peak curve. 

Remark 3.1. Usually, / 12 (ei) shouldbe small (f\ 2 {e\) = 
for f\ being linear). Thus, the sign ofG(e\,e 2 ) is primarily 
determined by the sign of f 22 (e 2 ). So, one may say that a 
necessary condition for the snap-through (i.e., there are two 
turning points in the post-peak curve) is the change of the 
convexity of the constitutive curve of the softening part. 

Now, let us consider the differences between Case B and 
Case C. Recall that e* n and e* n (e* u > e* n ) are the two roots 
of equation (13). For Case B, we have 



(I -n)e* 12 + ne2 2 < (I -n)eiQ + ne 20 , where e* 22 = f 2 [fi{e* l2 ) 

(27) 
On the other hand, 

g(e* n ,e* 22 ;n) = nf[{e\ 2 ) + (1 - n)f' 2 (e* 22 ) =0. (28) 

From the above two equations, we can get (26) immediately. 
Also, from (26) and (28) one can immediately deduce (27). 
This completes the proof. 

Assumption (25) can be made even more complicated, in that 
case we may draw the fairly complicated post-peak curves 
in 0, which were obtained by numerical methods. It should 
be pointed out that inequality (26) is another requirement 
among f\, f 2 and n. For given f\ and f 2 , it provides an- 
other bound (say n 2 ) for n, since e\\ and e| 2 are related to n 
as equation (28) shows. 

For f\ being linear in the post-peak region, it is easy to 
show that (26) becomes 



Meh) > 



/2O22) ~ /2(£20) 



'22 



e 2 o 



(29) 



(l-nK2+n/2 _1 [/i(^2)]>5o 



(23) 



This implies that for the constitutive relation O = f 2 (e 2 ), the 
secant line joining the point e* 21 and the peak e 2 ?> should be 
steeper than the tangent line at e 22 (see Figure 4). Combined 
with Corollary 3.1, we have the following corollary: 

Corollary 3.2. For fi being a linear function, if the con- 
vexity of f 2 changes once, then a necessary and sufficient 
condition for the occurrence of Case C is n G [«i,«o], and 
inequality (29) holds. 




the secant line 
the tangent line 



e2o e22 

Fig. 4. Diagrammatic representation of inequality (29) 



4 Illustrative Examples 

In this section, we give two examples to illustrate the 
theoretical results obtained in Section 3. The following two 
examples can be referred to as two different physical pro- 
cesses. One can easily check that the functions in the fol- 
lowing examples satisfy the conditions we have proposed, in 
particular, (25). 

Example 1: Consider the following constitutive rela- 
tions: 



971 



/n(ei)=csc — sin(507Cei), 



9tc 9tc 9tc 

fv(ei) = (507ccsc — )ei esc \- 1, 



3.3 Case D: Snap-back 

In Case D, there is a snap-back in the structural re- 
sponse. Here, we say that snap-back occurs when the slop of 
the force-displacement curve becomes positive and remains 
positive in the post-peak response. Obviously, in this case 
there is only one turning point (see Figure 2). The following 
theorem provides a critical n for the occurrence of Case D. 

Theorem 3.3. If G(ei,e2) changes sign once (cf(25)) or 
does not change sign for e\ £ [en,eio]> then for Case D to 
occur, a necessary and sufficient condition is n <n\. 

Proof. First, suppose that G{e\,e%) changes sign once. Then 
m(ei,e2) has the characteristics shown in Figure 3. It is ob- 
vious that a necessary and sufficient condition for the occur- 
rence of Case D is that there exists only one root for equa- 
tion (12). While from Figure 3, it is easy to see that a nec- 
essary and sufficient condition is n <n\. Second, suppose 
that G(e\,e2) does not change sign. Since dm/de2 > (as 
dm/de2 > at ei — e2o), we have 



«o = max m(ei,e2) =m(en,e2i)=ni. (30) 

e 2 G[e20,e2l] 



Obviously, a necessary and sufficient condition for n = 
m(ei,e2) to have one and only one root (i.e., g(e\,e2',n) has 
one and only one zero) is n < n\. Thus we complete the 
proof. 

Remark 3.2. In this section, we derive some requirements 
on the constitutive functions, together with three critical val- 
ues (no, n\, n%) of the size parameter. Providing the constitu- 
tive requirements are met, the structural response may have 
different behaviors for n in different intervals according to 
the above critical values. Thus, the results also show the 
important size effects. 



h\{e 2 ) =sin(507te 2 ), 



/22(e2)=65625[-(e 2 -0.03) 3 -0.0004(e 2 -0.03)]+0.65. 



Here Go = 1, and we take a* = 0.301. The details are shown 
in Figure 5(a). 

We have taken /12 being a linear function, which rep- 
resents the physical situation that at the peak region A of 
the structure has entered the plastic state. We find the crit- 
ical values of n based on the theoretical analysis in Sec- 
tion 3: « = 0.142, n\ = 0.0158, n 2 = 0.110 (a bound for 
n found from inequality (29)). Specifically, Case A occurs 
if n> 0.142; Case B occurs if 0.110 < n < 0.142; Case C 
occurs if 0.0158 <n< 0.110; Case D occurs if n < 0.0158. 
Accordingly, by taking n to be in different intervals, we get 
the four cases as we have predicted in Section 3. They are 
shown in Figure 5(b). 

To reflect the size effect on the localization zone, curves 
of the width of the localization zone in the current config- 
uration versus the total elongation are shown in Figure 6. 
This width is denoted by d, whose expression is given by 
d = n{\ +e2). Here, for the purpose of clearness, we have 
used different scales for different curves. It can be seen that 
this width increases slowly in the pre-peak region and in- 
creases rapidly in the post-peak region. For the stable soft- 
ening case (n = 0.167), there is only one value of d for a 
given e. For the snap-back case (« = 0.0156), there are two 
values of d for a given e. Also, d increases very fast, as e 
decreases in the post-peak region. For the two snap-through 
cases (n = 0.1 and n = 0.125), there are three values of d for 
e in some intervals. Thus, d may jump from a small value 
to a large value for e in these intervals, i.e., the localization 
zone may suddenly widen. Thus, the size parameter n has an 
important influence on the localization zone. 
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Fig. 5. (a) The constitutive curves of Example 1 ; (b)The engineering 
stress-strain curves for different n in Example 1 . 
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Fig. 7. (a) The constitutive curves of Example 2; (b) The engineer- 
ing stress-strain curves for different n in Example 2. 
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f2i(e 2 ) = -— ( e2 -0.009) 2 +l, 



fa{ei) =65625 [-(e 2 -0.029) 3 - 



0.0004(e 2 -0.029)]+0.65. 
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Fig. 6. The d — e curves for different n in Example 1 



Here Go = 1, and we take a* = 0.333. The details are 
shown in Figure 7(a). Critical values of n are: «o = 0.174, 
n\ = 0.0614, n 2 = 0. 136 (a bound for n found from inequal- 
ity (26)). The intervals for different cases are: Case A occurs 
if n > 0.174; Case B occurs if 0.136 < n < 0.174; Case C 
occurs if 0.0614 <n< 0.136; Case D occurs if n < 0.0614. 
The curves for n taking four values in these four different 
intervals are shown in Figure 7(b), which agree with our the- 
oretical predictions in Section 3. Curves of the width of the 
localization zone in the current configuration versus the total 
elongation are also shown (see Figure 8). Once again, from 
these curves, one can see the important influence of the size 
parameter n on the localization zone. 



Example 2: In this example, region A is assumed to be 
in nonlinear elasticity (loading or unloading), so the consti- 
tutive functions f\\ and fyi are the same. The constitutive 
relations are listed below. 



/n(ei)=/i2(ei) = -11000(ei-0.01) 2 +l.l, 



5 Concluding Remarks and Future Tasks 

An analytical study is performed on the post-peak struc- 
tural response of strain-softening with localization. In a 
general nonlinear setting, after taking standard mathematical 
analysis to the parametric equations, we manage to handle 
the nonlinear and size effects. Qualitative requirements on 
the constitutive functions and quantitative requirements on 
the size effect are derived, especially for the snap-through 
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phenomenon. The results are consistent with earlier experi- 
mental and computational results. It seems that the four cases 
studied analytically here are quite representative. The theo- 
retical results may be of value for the verification of compu- 
tational algorithms and can shed some light on the mecha- 
nisms of instabilities associated with strain-softening. Espe- 
cially, we have shown that the convexity change is a neces- 
sary condition for the snap-through phenomenon. As soften- 
ing with localization is important for understanding the fail- 
ure evolution in structures, future work will focus on consid- 
ering structures with different configurations. 
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